UNIVERSITY 
OF  WISCONSIN 


CENTER  FOR  THE 
MATHEMATICAL 
SCIENCES 


Center  for  the  Mathematical  Sciences 
University  of  Wisconsin — Madison 
1308  W.  Dayton  Street 
Madison,  Wisconsin  53715-1149 


February  1996 
(Received  February  9,  1996) 


19960227  113 

Approved  for  public  release 
Distribution  unlimited 


Sponsored  by 


National  Science  Foundation 
4201  Wilson  Blvd. 

Arlington,  VA  22230 


U.  S.  Army  Research  Office 
P.O.  Box  12211 

Research  Triangle  Park,  NC  27709 


UNIVERSITY  OF  WISCONSIN-MADISON 
CENTER  FOR  THE  MATHEMATICAL  SCIENCES 

Compactly  supported  tight  affine  spline  frames  in  L2(IR*^) 


Amos  Ron 


Zuowei  Shen 


Computer  Science  Department 
University  of  Wisconsin-Madison 
1210  West  Dayton  Street 
Madison,  Wisconsin  53706,  USA 
ainos@cs.wisc.edu 


Department  of  Mathematics 
National  University  of  Singapore 
10  Kent  Ridge  Crescent 
Singapore  119260 
matzuows@leonis.nus . sg 


Technical  Summary  Report  #  96-0  5 

February  1996 

ABSTRACT 

The  theory  of  [RS2]  is  applied  to  yield  compactly  supported  tight  affine  frames 
(wavelets)  in  from  box  splines.  The  wavelets  obtained  are  smooth  piecewise- 

polynomials  on  a  simple  mesh;  furthermore,  they  exhibit  a  wealth  of  symmetries,  and  have 
a  relatively  small  support.  The  number  of  “mother  wavelets”,  however,  increases  with  the 
increase  of  the  required  smoothness. 

Two  bivariate  constructions,  of  potential  practical  value,  are  highlighted.  In  both,  the 
wavelets  are  derived  from  four-direction  mesh  box  splines  that  are  refinable  with  respect 

to  the  dilation  matrix 
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Compactly  supported  tight  affine  spline  frames  in  L2(1R'^) 
Amos  Ron  and  Zuowei  Shen 


1.  Introduction 

Given  a  finite  set  ^  C  L2(1R‘^),  and  a  dilation  matrix  s,  the  affine  system  generated 
by  is  defined  as  the  collection 

(1.1)  A  :  V'  e  A;  G  K,  a  G  2Z‘^}, 

where 

E^:f^f{-  +  a) 

is  the  shift  operator  and 

is  the  dilation  operator.  A  system  X  C  L2(IR‘^)  is  a  fundamental  tight  frame  with 
frame  bound  1  if  the  map 


r*  :  L2(IR‘')  ^  {{f,x)),^x 

is  unitary  (but  not  necessarily  onto:  a  tight  frame  with  frame  bound  1  whose  corresponding 
T*  is  onto  is  necessarily  orthonormal).  In  what  follows,  all  systems  that  we  treat  are  affine 
and  fundamental,  and  all  tight  frames  that  are  considered  have  frame  bound  1,  hence  “a 
tight  frame”  should  always  be  understood  as  “a  fundamental  tight  affine  frame  with  frame 
bound  1” .  A  tight  frame  can  be  used  for  the  atomic  decompositions  of  functions  exactly 
in  the  same  way  orthonormal  bases  are  used,  i.e., 

T*:f^  T*f 

transform  /  into  discrete  information,  and  its  adjoint 

T  :  ^2(A)  ->  L2  :  c  i->  ^  c{x)x 

xex 

can  then  be  used  to  recover  /  from  its  discrete  transformation  (i.e.,  TT*f  =  /).  In  order 
for  T*  to  exhibit  good  space-frequency  localization,  the  functions  should  be  local  in  the 
space  domain  (ideally  compactly  supported),  should  be  smooth  (which  leads  to  good  decay 
in  the  frequency  domain),  and  should  provide  positive  approximation  order.  We  refer  to 
[Dl,2],  [HW]  and  [RSI, 2]  for  further  discussions  on  frames  and  tight  frames. 

While  tight  frames  should  ideally  be  generated  by  few  smooth  compactly  supported 
functions  with  simple  structure,  the  only  examples  of  compactly  supported  multivariate 
tight  frames  in  the  literature  that  we  are  aware  of  are  bivariate  orthonormal  bases  that 
were  derived  in  [CD]  from  the  univariate  Daubechies’  wavelets,  [D3],  as  well  as  tensor 
products  of  Daubechies’  wavelets.  We  briefly  review  those  constructions  now. 
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Let  L2  :=  L2(]R^).  Let  h  be  the  refinement  mask  of  a  Daubechies’  scaling  function  (j), 
[D3]  (whose  shifts  are  known  to  be  orthonormal),  and  let  tp  be  the  corresponding  wavelet. 
Define  the  bivariate  mask 

t{ui,U2)  :=  h{uJi). 

Given  a  dilation  matrix  s,  a  bivariate  scaling  function  (or  distribution)  $  can  then  be 
defined  by 

00 

$  :  w  J][ 

i=i 

Cohen  and  Daubechies  employed  in  [CD]  two  different  dilation  matrices: 


and  obtained  therefore  two  different  scaling  functions,  say  $  and  $1,  respectively.  They 
observed  essential  differences  between  the  two  so-obtained  functions.  The  matrix  s  sat¬ 
isfies  5’^  =  27,  and  therefore  the  resulted  $  is  clearly  the  tensor  product  (referred  to  as 
“separable”  in  [CD]) 

^{x)  -  (I){x2)(p{xi  -  3:2), 

and  therefore  the  shifts  of  $  are  necessarily  orthonormal.  Furthermore,  the  standard 
wavelet  construction  then  yields  the  wavelet 

’^(r;)  =  ip{x2)(l)(,xi  -  X2). 

While  the  refinable  function  $  is  separable,  the  refinable  function  and  the  wavelet 
constructed  from  it  are  not  separable  in  any  sense.  Nonetheless,  it  is  proved  in  [CD]  that 
the  shifts  of  $1  are  orthonormal.  It  is  further  proved  in  [CD]  that,  unfortunately,  $1 
cannot  be  regardless  of  the  order  of  the  univariate  scaling  function  which  is  used. 
Therefore,  it  is  necessary  to  develop  some  other  algorithms  to  construct  nonseparable 
compactly  supported  tight  frames  with  high  smoothness. 

It  will  not  be  entirely  correct  to  say  that  the  [CD]  constructions  and  tensor  product 
constructions  comprise  all  known  multivariate  affine  tight  frames:  according  to  [CS] ,  tight 
frames  can  be  constructed  by  appending  to  some  of  their  translates  (=:oversampling); 
that,  of  course,  only  increases  the  number  of  elements  in  each  of  the  above  constructions, 
while  preserving  any  deficiences  (such  as  lack  of  symmetries  and  low  smoothness  in  the 
non-separable  case,  and  parallelogram  supports  in  the  separable  case)  that  the  orthonormal 
system  may  have  had. 

The  theory  established  in  [RS2],  however,  makes  the  construction  of  useful  simple 
nonseparable  compactly  supported  tight  frames  with  high  smoothness  and  a  variety  of 
symmetries  an  easy  task.  In  fact,  one  can  essentially  construct  tight  frames  with  the  aid  of 
the  shifts  of  any  refinable  function,  and  therefore  it  is  possible  to  impose  a  simple  structure 
on  the  wavelets  by  selecting  a  refinable  function  with  such  desired  structure.  Indeed,  we 
constructed  in  [RS2],  for  every  positive  integer  m,  a  tight  frame  for  L2(IR)  that  is  generated 
by  m  wavelets,  each  of  which  is  a  spline  of  degree  m  -  1,  support  [0,  m],  and  smoothness 
(^m-2  Further,  all  knots  of  the  spline-wavelets  are  half-integers,  and  each  spline  is  either 
symmetric  or  anti-symmetric.  (The  case  m  =  4  of  this  construction  is  discussed  in  §4). 
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Since  the  variety  of  possible  constructions  of  multivariate  (and  univariate)  tight  frames 
based  on  the  [RS2]  theory  is  unlimited,  we  carefully  selected  for  the  present  article  those 
constructions  which,  in  our  opinion,  may  be  used  in  practical  applications.  We  kept  in 
mind  that  different  applications  may  require  different  properties  from  the  wavelet  system; 
for  example,  in  data  compression  applications  the  number  of  different  wavelets  used  (which 
accounts  to  the  oversampling  rate)  should  be  minimized,  while  in  finite  element  applica¬ 
tions  many  elements  with  small  supports  may  be  preferred.  Our  two  favorite  bivariate 
constructions  are  detailed  in  §2  and  §3:  in  §2,  the  system  is  generated  by  many  highly 
symmetric  wavelets  of  small  support  (which,  actually,  are  not  so  “many”  for  practical 
smoothness  requirements);  in  §3,  the  system  is  generated  by  fewer  larger  elements.  In  each 
construction,  the  wavelets  are  splines,  i.e.,  smooth  piecewise-polynomials,  always  of  com¬ 
pact  support,  with  various  symmetries.  Moreover,  the  relevant  grid  is  the  four-direction 
(i.e.,  quincunx,  see  below)  mesh,  hence  all  the  wavelets  have  a  “round”  octagonal  support. 
On  the  the  other  hand,  similarly  to  our  univariate  construction,  the  number  of  elements 
used,  as  well  as  the  volume  of  the  support  of  each,  increase  together  with  the  increase  in 
the  required  smoothness.  The  ease  in  constructing  tight  spline  frames  is  not  limited  to  two 
dimensions:  in  §4  a  general  inductive  algorithm  for  constructing  multidimensional  tight 
frames  is  provided.  The  algorithm  works  particularly  well  with  box  splines,  in  fact,  with 
every  box  spline.  We  call  the  box  spline  wavelets  obtained  by  this  algorithm  boxlets. 

We  have  chosen  to  carry  the  constructions  in  §2,3  with  respect  to  the  dilation  matrix 
s  (cf.  1.2).  It  is  possible,  though,  to  carry  these  same  constructions  with  respect  to  si 
with  only  one  limitation:  all  four  directions  in  the  definition  of  the  box  spline  must  appear 
with  the  same  multiplicity.  In  any  event,  the  use  of  S\  instead  of  s  does  not  yield  different 
systems,  and  the  reason  is  very  simple:  for  a  four-direction  mesh,  not  only  that  sTL^  — 
siTZ^,  but  also,  in  case  the  box  spline  (f)  has  equi- multiplicities,  0(s-)  =  ^(si-). 


2.  Bivariate  tight  frames  generated  by  —  1  ^  spline  wavelets 

In  fact,  the  construction  here  is  a  bit  more  general  than  the  title  indicates:  it  includes 
tight  frames  generated  by  —  1  wavelets  of  smoothness  Of  course,  we  assume 

m  >  1. 

Let  (p  be  the  box  spline 


(2.1) 


= n 


■  U) 


where 


(^15^2,  ^3»^4) 


5 


and  with  mi  =  m3,  and  m2  =  rn^.  It  is  well-known  (cf.  [BHR])  that  the  box  spline  0 
satisfies  the  following  properties: 

(a)  It  is  a  piecewise-polynomial  of  local  degree  2(mi  -l-m2  —  1),  on  the  four-direction  mesh 
(that  is  obtained  by  adding  the  diagonals  to  each  square  of  integer  vertices). 
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(b)  It  is  globally  ^  with  p  :=  min{2mi  +  m2,  2m2  +  mi};  (in  fact,  its  (p  2)-order 
derivatives  are  all  Lip^,  hence  its  Holder  continuity  is  p  -  1).  Further,  it  provides 
approximation  order  p. 

(c)  It  is  supported  in  the  octagon 

4 

^3^3  '■  0  <  J  =  1,  •  •  •  5  4}. 

:/=i 

Further,  it  is  essentially  positive  on  its  support. 

In  the  sequel,  we  will  introduce  a  variety  of  octagonal  domains  similar  to  the  support 
of  0  above.  For  reasons  of  efficiency,  we  therefore  denote 

4 

[oi,  0,2, 03,  ^4]  •=  •  0  ^  —  %■)  j  —  ^1-  ••  ■> 4}- 

j=i 

Note  that  each  [a,  b,  c,  d]  is  an  octagon  whose  area  is 

area([o,  6,  c.,d])  =  ab  +  ac  +  ad  +  be  +  bd  +  2cd. 


We  prefer  the  above  box  spline  over  other  variants  because  it  is  refinable  not  only 
with  respect  to  the  dilation  matrix  27,  but  also  with  respect  to  the  dilation  matrix 


(2.2) 


This  well-known  fact  was  rarely  exploited  before  because  the  shifts  of  the  4-direction  box 
spline  do  not  form  a  Riesz  basis  (or  a  frame).  Here,  that  is  no  more  an  obstacle,  since,  as 
we  already  mentioned,  the  [RS2]  theory  does  not  impose  any  a  priori  assumption  on  the 
refinable  function  (j). 

A  straightforward  computation  shows  that 

^{S-)  =  To(j) 


with  the  function 


Toico) 


/ 1  + 

V  2 


mi 


1  -F 

2 


m2 

=  0{u)  cos”*‘  {u3i/2)  cos^^  {U2/2), 


with 

+Tn2ii;2)/2^ 

We  note  that  4>  is  refinable,  which  means  that  the  mask  tq  is  27r-periodic  (and  hence 
(f){s~^-)  can  be  written  as  a  linear  combination  of  the  integer  translates  of  (j)). 

Our  construction  invokes  the  following  theorem,  which  is  a  special  case  of  Corollary 
6.7  of  [RS2]. 
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Theorem  2.3.  Let  (f)  be  any  box  spline  that  is  rehnable  with  respect  to  the  dilation  matrix 
s  of  (2.2),  and  has  a  (2n-periodic)  refinement  mask  tq.  Let  n  be  a  positive  integer,  and  let 
(Tj)^_i  be  n  27r-periodic  essentially  bounded  measurable  functions.  Assume  that,  for  a.e. 
u),  and  for  v  :=  (tt,  tt), 

Ekp  =  i.  j2r,wy=(s. 

3=0  3=0 

Then  the  wavelets  ^  :=  defined  by 


j  =  l,...,n, 

generate  a  tight  frame  (i.e.,  a  fundamental  affine  tight  frame  with  frame  bound  1)  for  L2. 

It  should  be  understood  that  the  dilation  matrix  s  that  is  involved  in  the  refinement 
equation  is  the  same  dilation  matrix  s  that  is  used  to  generate  the  affine  system  from  '5'. 

To  construct  now  the  wavelets  with  the  aid  of  the  above  box  spline  (j>,  we  first  define 
the  following  univariate  47r-periodic  functions; 


(2.4)  yj{t)  :=  yf^\t)  :=  y  cos’"  ^  t/2  sin^  t/2,  0  <  j  <  m. 

(Warning;  m  in  is  an  index,  not  a  power!)  Note  that,  for  any  fixed  m, 

m  m 

(2.5)  =  ^Z/i2/j(-  +  7r)  =  0. 

3=0  j=o 

We  next  define  the  tensor  product  bidimensional  mask  system 

(r„  ;  n  G  W  ;=  {0,l,...,mi}  X  {0,l,...,m2}) 

with 

Each  Tn  is  27r-periodic. 

We  then  observe  that  tq  ;=  ro,o  in  the  above  tensor  system  is  the  refinement  mask  of 
(j).  Moreover,  the  convex  hull  of  the  spectrum  of  each  r^,  n  G  N,  is  independent  of  n  (i.e., 
in  down-to-earth  language;  the  27r-integers  that  lie  in  the  convex  hull  of  the  frequencies  of 
the  exponentials  whose  linear  combination  form  r„,  are  the  same  for  all  n).  This  means 
that  the  (mi  -|-  l)(m2  -f- 1)  —  1  functions,  defined  by 

(2.6)  ;=  {tpn  ■  ^n(s-)  :=  n  6  iV\0}, 

are  all  supported  in  the  support  of  (f),  i.e.,  in  the  octagon  [mi,  m2,  mi,  m2].  The  functions 
4/  are  nonseparable  piecewise-polynomials,  with  the  same  smoothness  as  the  box  spline  (j) 
used. 
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The  tensor  product  structure  of  the  masks  (together  with  the  fact  that  |0|  —  1)  clearly 
implies  (cf.  (2.5))  that 


mi 

neN  j-0  j=0 

and  by  the  same  token,  for  u  =  (tt,  tt)  (with  c  :=  9  E''6) 

TUi  _ _ _  - - 

E  T„(w)W^)  =  c(E!/j“'>(a.i)!,!™-l(a,,  +ir))(E!/l™'(‘^2)!/‘"*“'(w2  +  >r))  =  0-0  =  0. 

n£N  i=0  i=0 

Therefore,  Theorem  2.3  implies  the  following: 

Theorem  2.7.  The  wavelets  ^  constructed  in  (2.6)  generate  a  tight  affine  frame  with 
frame  bound  1. 

Remark.  Each  one  of  these  wavelets  is  supported  in  the  support  of  <f),  and  is  a  C''’ 
p  :=  min{2mi  +m2, 2m2  +  mi},  piecewise-polynomial  of  local  degree  2(mi  +  m2  - 1),  with 
respect  to  the  mesh  s~^M,  with  M  the  standard  four-direction  mesh.  The  mesh  s  M 
consists  of  all  lines  of  the  form 


X\-\-X2=  j,  Xi  -  X2=  j,  Xi  —  i/2,  X2  ~  i/2, 
where  j  varies  over  71,  and  {xi,X2)  is  the  generic  point  in  the  space  domain. 


Example.  We  take  mi  =  m2  =  1-  In  this  case  0  is  the  well-known  piecewise-quadratic 
Zwart  element,  supported  in  the  octagon  [1,1,1,!].  The  above  construction  yields  three 
piecewise-quadratic  wavelets  each  supported  in  [1, 1, 1, 1],  and  each  exhibiting  various 
symmetries  (cf.  Figures  2.1  and  2.2).  The  support  of  each  wavelet  is  of  area  7,  this  small 
support  comes  despite  of  fact  that  we  oversample  by  a  factor  of  3. 

It  is  hard  to  compare  the  above  construction  to  literature  counterparts,  since,  as  we 
mentioned  before,  the  latter  hardly  exist.  Here  are  two  possible  comparisons.  In  the 
first,  we  take  the  separable  construction  of  [CD],  using  Daubechies  scaling  function  that 
provides  the  same  approximation  order  (i.e.,  3)  as  the  Zwart  element.  The  result  is  a 
single  wavelet  supported  in  a  parallelogram  of  area  25,  and  whose  shifts  are  orthonormal. 
The  smoothness  of  this  wavelet  is  about  half  of  the  above  tight  frame  (in  terms  of  Hdlder 
exponents:  2  vs.  1.1;  cf.  [D2:  pp.  232-239]),  and  it  lacks  symmetries. 

Another  comparison  may  be  with  “plain”  tensor  product  of  the  same  Daubechies 
scaling  function.  In  this  case,  we  obtain  3  mother  wavelets  each  with  a  square  support 
of  area  25,  and  with  the  same  smoothness  and  (lack  of)  symmetries.  The  shifts  of  these 
wavelets  now  fill  an  entire  dyadic  level.  To  fill  in  an  entire  dyadic  level  using  integer  shifts 
of  our  elements,  we  need  the  three  elements  of  support  area  of  7  each,  and  6  elements  of 
half-size  support,  with  total  area  of  supports  42.  Id 
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Figure  2.1.  The  (7^  piece wise-quadr atic  wavelet '0o,i-  The  wavelet 
V'l.o  is  obtained  by  rotation. 


Figure  2.2.  The  piecewise-quadratic  wavelet  •0i,i.  Note  the 
octagonal  support 

Example.  We  take  mi  =  1,  m2  =  2.  In  this  case  (j)  is  globally  (and  is  on  half  of 
the  mesh  lines),  piecewise-quartic.  The  construction  now  yields  5  wavelets,  all  supported 
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in  the  octagon  [1,2, 1, 2]  whose  area  is  15.  They  all  have  the  same  (global)  smoothness  and 
local  degree  as  the  box  spline. 

Example  2.8.  We  take  mi  =  m2  =  2.  Then  we  obtain  8  elements  supported  each  in  the 
octagon  [2, 2, 2, 2]  whose  vertices  (up  to  a  (1,  l)-shift)  are  (±3,  ±1)  U  (±1,  ±3)  and  whose 
area  is  28.  Each  spline  wavelet  is  (with  Holder  exponent  5).  The  local  polynomial 
degree  is  6. 


3.  Tight  frames  generated  by  2m  ^  spline  wavelets 


Here,  m  assumes  any  positive  integer  value.  Again,  the  construction  is  more  general 
than  the  title  indicates:  it  includes  the  construction  of  tight  frames  generated  by  2m  +  1 
(jZm-i  vvavelets. 

The  wavelets  constructed  in  this  section  are  selected  from  the  same  box  spline  spaces 
used  in  the  previous  section.  Given  a  four-direction  box  spline  with  corresponding  mul¬ 
tiplicities  (r?2i,m2,mi,m2),  the  construction  selects  appropriate  mother  wavelets  as  the 
s-dilate  of  certain  functions  in  the  span  of  the  shifts  of  the  box  spline.  This  means  that 
the  wavelets  here  are  comparable  to  the  wavelets  of  the  previous  section  in  terms  of  local 
polynomial  degrees,  underlying  meshes,  and  smoothness. 

However,  the  number  of  wavelets  associated  with  the  previous  construction  grows 
quadraticly  with  the  required  smoothness.  Though  the  three  examples  that  followed  show 
that  for  practical  smoothness  requirements  the  number  of  wavelets  is  “within  reason” ,  it 
is  possible  to  construct  tight  spline  frames  whose  number  of  generators  grow  only  linearly 
with  the  required  smoothness.  The  associated  wavelets,  on  the  other  hand,  have  larger, 
somewhat  less  symmetric,  support. 

It  is  convenient  to  carry  out  the  construction  here  not  with  the  aid  of  the  box  spline 
(f)  of  (2.1),  but,  rather,  with  the  following  “averaged”  spline  cp: 


(p{w)  := 


1  + 

2 


mi 

^(cu). 


In  standard  box  spline  terminology,  is  a  box  spline  with  direction  set  (2^i,^2^3)^4)  (cf- 
(2.1))  with  corresponding  multiplicities  (mi,  m2,  mi,  m2).  Also,  direct  computation  yields 
that  (f  satisfies 


(p{s(jj)  = 


1  +  /  I  + 


-IU)2 


m2 


^(w) 


-d{u})  cos"^^  ((wi  +  a;2)/2)  cos”^^  (^>2/2)  ^(u>) 


where 


_ g— z(mi  ,mi+m2)-tt;/2 


Recalling  now  the  definition  of  in  (2.4),  we  define  mi-f  m2+l  27r-periodic  functions 
as  follows: 

Tj{u))  -f  0^2)  2/0^^  (^2)5  j  =  0, . . .  ,mi. 
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and 

Tm^+j{u})  :=  9{uj)yf^^\iV2),  i  =  1, . . .  ,m2, 

where 

e{uj)  := 

We  note  that  the  functions  are  periodic  with  respect  to  shifting 

by  (7r,7r),  hence  are  periodic  with  respect  to  the  lattice  s~^2'k7L^.  From  that,  together 
with  the  fact,  (2.5),  that 

m\ 

Ei»f‘'P  =  ‘. 

3=0 

one  immediately  concludes  that  for  £  €  {0, 1},  and  =  (tt,  tt), 

771 1  - 

=  6{-£'k)  +  ^Tt). 

3=0 

This  leads  further  to  the  conclusions  that,  firstly, 

7711+7712  7712 

E  l’■.l^  =  El!'^'P  =  l. 

3=0  3=0 

and,  secondly,  with  v  :=  (7r,7r), 

7711+7712  7712  - 

rj{u3)E''Tj{u3)  =  9{-'K)^^yY^\uJ2)yY^\^‘^  +  tt)  =  0. 

3=0  3=0 

In  summary,  the  mask  vector 

\‘mi+m2 

^3)3=0 

satisfies  the  conditions  of  Theorem  2.3. 

Consequently,  since  tq  in  the  above  construction  is  the  refinment  mask  of  </?,  we  con¬ 
clude  from  Theorem  2.3  that  the  mi  -|-  m2  wavelets  defined  by 

ijjj{s-)  Tj(p,  j  —  1, . . .  ,mi  +  m2, 

generate  a  tight  frame. 

Corollary  3.1.  Let  ip  be  the  box  spline  whose  direction  set  is  (2^i,^2)^3)^4)  with  corre¬ 
sponding  jjiultiplicities  (mi,  m2,  mi,  m2): 

- 2 - J 
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with  (t>  the  box  spline  of  (2.1).  Then,  with  the  above  constructed  masks,  the 

wavelets  ^  :=  defined  by  $j{s-)  :=  Tjip,  generate  a  tight  frame  (with  frame 

bound  1)  for  L2-  These  wavelets  share  the  same  smoothness,  local  degree,  and  mesh  with 
the  wavelets  of  Theorem  2. 7. 


Discussion.  The  box  spline  ip  in  the  above  corollary  is  supported  in  the  octagon 
[2mi, m2, mi,  m2]  whose  area  is  2m\  +  7m\m2  +  rn^.  As  to  the  wavelets,  the  “large” 
wavelets,  i.e.,  the  first  mi  elements,  have  the  same  support  as  (p.  The  other  m2  elements 
have  the  octagonal  support  [mi,  m2,  mi,m2\,  whose  area  is  m\  +  5mim2  + 

Remark.  The  general  construction  detailed  in  the  next  section  shows  that  there  are  many 
possible  modifications  of  the  above  construction. 

Example  3.2.  We  consider  the  case  mi  =  m2  =  1-  Figure  3.1  shows  the  box  spline  ip.  It 
is  a  piecewise-quadratic  that  provides  approximation  order  3,  and  it  is  supported  in  the 
octagon  [2, 1, 1, 1],  with  area  10.  Two  wavelets  are  constructed  here.  Direct  computation 
yields  that  the  larger  support  wavelet,  ipi  (cf.  Figure  3.2),  whose  support  is  identical  to 
that  of  ip,  satisfies 

(3.3)  ^1(0;)  =  tan(a;i/2)^(a;), 

while  the  smaller  support  ■02  (cf.  Figure  3.3)  has  the  form 


*  = - 

The  corresponding  supports  are  [2, 1, 1, 1]  and  [1, 1, 1, 1]  with  corresponding  areas  10  and 


7. 


□ 


Example.  We  consider  mi  =  1,  m2  =  2.  The  box  spline  ip  is  then  piecewise-quartic 
and  provides  approximation  order  4.  Its  support  is  the  octagon  [2,2, 1,2]  whose  area  is 
20.  The  mother  wavelet  set  consists  of  three  elements,  two  “small”  and  one  “large” .  The 
large  element  0i  is  defined  exactly  as  in  (3.3)  (only  that  ip  has  been  changed),  and  has  the 
support  of  (j).  The  other  wavelets  are 

^  1  —  ”■“'2 1  ^  u)\  +  Ct<2  ^1  ”  ^2  >. 


and 


^3  = 


-1  +  2e- 


^  +  ^2  ijOi  —  UJ2 


-ip( 


), 


Both  supported  in  the  same  octagon  [1, 2, 1,  2],  whose  area  is  15. 


□ 


Example.  The  case  mi  =  m2  =  2  yields  four  wavelets,  two  supported  in  [4, 2, 2, 2],  and 
the  other  two  supported  in  [2, 2,  2,  2].  The  area  of  these  octagons  is  40  and  28  respectively, 
hence  their  total  area  is  136.  The  smoothness  of  the  box  spline,  as  well  as  its  local  degree 
is  the  same  as  in  Example  2.8. 
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Figure  3.1.  The  piecewise-quadratic  box  spline  that  is  used 
in  Example  3.2. 


Figure  3.2.  The  piecewise-quadratic  wavelet  with  larger  sup¬ 
port  ipi. 
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0.2^ 


Figure  3.3.  The  piecewise-quadratic  wavelet  with  smaller  sup¬ 
port  'tp2- 

4.  A  general  inductive  algorithm  for  constructing  tight  frames 

The  method  described  in  the  last  section  for  the  construction  of  box  spline  wavelets 
on  a  four-direction  mesh  can  be  significantly  generalized.  Since  we  are  unable  to  predict, 
at  the  time  when  this  article  is  written,  what  specific  variants  may  be  implemented  in 
practice,  we  decided  to  simply  outline  the  highest  level  of  generalization  that  we  are  able 
to  observe. 

The  setup  here  is  as  follows:  we  hold  a  dilation  matrix  s  (i.e.,  an  integer  matrix  whose 
inverse  is  contractive)  and  two  d- variate  functions  (or,  more  generally,  distributions),  say 
</>i,  (f)2,  that  are  refinable  with  respect  to  s.  The  corresponding  masks,  that  are  assumed 
to  be  27r-periodic  (that  is  a  part  of  the  definition  of  refinability),  are  also  assumed  to  be 
bounded;  in  practical  situations  the  masks  are  trigonometric  polynomials,  hence  certainly 
bounded.  The  basic  assumption  is  that  we  already  know  how  to  derive  a  tight  frame  from 
the  scaling  function  <^i,  and  we  would  like  to  use  that  known  frame  in  order  to  obtain  a 
new,  improved  (in  terms  of  smoothness,  for  example)  tight  frame.  Specifically,  we  would 
like  to  extract  the  new  frame  from  the  (necessarily  refinable)  convolution 


(j)  •■=  (f>l*  <t>2- 


This  idea  (of  convolving  the  given  scaling  function  with  a  suitably  chosen  distribution) 
has  been  used  in  the  spline  and  wavelet  theory  many  times.  For  example,  Daubechies 
obtained  her  univariate  refinable  functions  whose  shifts  are  orthonormal  by  convolving  a 
univariate  B-spline  (f>\  (whose  shifts  are  stable  but  not  orthonormal),  with  a  suitably  chosen 
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distribution  ^2-  In  her  construction,  the  new  (j)  provides  the  same  approximation  order  as 
the  B-spline  but  is  significantly  less  smooth;  (the  smoothness  of  the  scaling  function 
and  the  orthonormality  of  its  shifts  are  then  transferred  to  the  constructed  wavelets.) 
The  fact  that  one  cannot  obtain  by  convolution  the  orthonormality  of  the  shifts,  while 
simultaneously  improving,  or  at  least  retaining,  the  smoothness  of  the  original  <^i,  is  one 
of  main  reasons  the  construction  of  the  affine  orthonormal  bases  is  fairly  involved. 

The  theory  of  [RS2]  allows  us  to  pay  less  attention  to  the  properties  of  the  scaling 
function:  “bad”  properties  of  the  latter  may  not  at  all  be  inherited  by  the  wavelets! 
Specifically,  the  smoothness  of  the  wavelets  and  the  tightness  of  the  frame  they  generate 
do  not  compete  any  more.  In  fact,  we  can  separate  the  construction  of  the  smooth  tight 
frame  into  two  steps:  in  the  first  a  basic  low-smoothness  frame  is  constructed,  and  then 
the  smoothness  of  the  frame  is  improved  without  hampering  its  tightness:  in  terms  of  0i, 
(f>2  above,  4>i  is  the  low-  (or  intermediate-)  smooth  scaling  function  that  yields  a  tight  affine 
system,  and  (j)2  is  the  convolutor  that  should  improve  the  smoothness. 

We  specify  now  the  exact  conditions  that  are  imposed  on  (f)i,  <j)2-  In  these  conditions, 
we  use  the  notation 

r  := 

For  example,  if  s  =  2/,  T  is  the  group  (0,  with  addition  modulo  27r.  Note  that  the 
order  of  F  is  |  dets|. 

(4.1)  Conditions  assumed  on  (pi,  (f)2- 

(a)  For  j  =  1,2,  (pj  is  a  refinable  distribution  with  bounded  mask  Tj. 

(b)  There  exists  a  collection  Ti  of  27r-periodic  bounded  functions  that  satisfy,  for  every 
i/  6  r,  and  a.e.  on  IR*^, 

tiE''ti  +  tE^'t  =  5i,. 

t€Ti 

(c)  There  exist  another  collection  of  27r-periodic  bounded  functions  T2,  that  satisfy 

\t2?  +  X! 

TeT2 


Discussion.  The  fact  that  we  assume  (p2  to  be  merely  a  distribution  seems  to  be  practically 
important  as  we  will  see  in  the  examples  presented  at  the  end  of  this  section.  On  the  other 
hand,  (pi,  in  all  examples  we  carry  in  mind,  is  a  function.  In  fact,  if  (pi  is  a  function. 
Corollary  6.7  of  [RS2]  “almost”  implies  that,  under  the  condition  (b)  of  the  conditions 
(4.1),  the  system  generated  by  the  wavelets  that  are  defined  via  their  Fourier  transforms 

by 

{'0(s*-)  :=  T01  :  r  e  Ti} 

is  a  tight  frame  for  L2(IR^).  The  only  missing  condition  is  a  very  mild  smoothness  re¬ 
quirement  of  (pi,  that  all  box  splines,  for  example,  satisfy  (cf.  (4.6)  of  [RS2]).  Note  that 
significantly  weaker  conditions  are  assumed  on  T2. 
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With  the  assumptions  of  (4.1)  in  hand,  we  will  be  able  to  construct  a  tight  frame 
based  not  on  the  convolution  ({>i  *  (j)2,  but,  rather,  on  the  convolution  with  larger  support 

|dets|  ■ 

Note  that  0  =  02(5*-)0i,  and  hence  that  0  is  refinable  with  mask 


:=  T2{s*-)  Ty. 


We  now  introduce  the  new  mask  collection 

T:=  TiU(riT2(s*-)), 

i.e.,  we  apply  s*  dilation  to  each  of  the  masks  in  T2,  multiply  each  by  ri,  and  then  append 
the  new  collection  to  Ti.  Clearly, 


#T  =  #Ti  +  #T2. 

Lemma  4.2.  Each  of  the  masks  in  T  above  is  2Tr-periodic  and  bounded.  Furthermore, 
for  every  n  G  T, 

rST 

Proof.  We  note  that  the  masks  in  T2(s*-)  are  not  only  27r!Z‘^-periodic,  but  also 
27rs*“^Z‘^-periodic.  Therefore,  for  every  €  F,  by  (c)  of  (4.1), 

Y.  (nr)E'^(riT)  =  {ryWTy)  (|r2(s*-)p  +  Y  kl")  = 

r€T2{s*)  tGT2(s*) 

This  implies  that,  by  (b)  of  (4.1),  that,  for  every  uGT, 

t^Wt^  +  Y  =  ryWTi  +  Y  ^ 

rGT  rgTi 

□ 

As  we  already  mentioned  above,  the  properties  proved  with  respect  to  T  in  the  above 
lemma  almost  imply  that  the  wavelet  system  ^  :=  (V’)v.)e^'  define  by 

{^(s*-)  :=  r0  :  r  G  T} 


generate  a  tight  frame;  we  only  a  need  a  mild  smoothness  condition  of  0.  Rather  than 
quoting  the  complicated  condition  (4.6)  of  [RS2],  we  assume  in  the  following  corollary 
smoothness  conditions  that  are  slightly  stronger  but  more  transparent. 
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Corollary  4.3.  In  the  above  notations,  if  the  rehnable  (f>  is  a  box  spline,  or  if  its  Fourier 
transform  decay  at  oo  like  0{\  •  for  5  >  0,  then  the  wavelets  ^  constructed  as 

above  generate  a  fundamental  affine  tight  frame  with  frame  bound  1  for  L2(1R‘^). 


(4.4)  Boxlets.  We  assume  that  the  dilation  matrix  s  satisfies  =  21,  for  some  positive 
integer  k.  In  this  case,  there  is  a  natural  way  for  choosing  the  smoothing  factor  (^2:  starting 
with  a  cycle  H  :=  of  s  (i.e.,  s^j  =  ^j+i,  j  =  I, ...  ,k  -  I,  hence  necessarily 

s^k  =  2^1),  we  define  ^2  as  the  box  spline  with  directions  B; 


= n 

^€3 


1  - 

i^-u) 


Then,  (/>2  is  refinable  with  mask 


'7*2  (ca)  = 


Since  we  tacitly  assume  B  to  be  integer,  T2  is  27r-periodic.  Selecting  T2  is  then  trivial:  T2 
is  taken  to  be  the  singleton 

T2  :={  2  7- 

The  satisfaction  of  (c)  in  (4.1)  is  then  automatic,  and,  assuming  that  Ti,  (f>i  are  already 
given,  and  that  the  refinement  mask  of  (pi  is  ri,  our  general  inductive  step  reads  here  as 
follows: 

(a)  The  new  refinable  function  is  (p  := 

(b)  The  refinement  mask  of  (p  is 

T^{u)  :=  Ti{oj) - - - =  Ti{u) - - - . 

(c)  The  new  wavelet  mask  set  T  is  obtained  by  appending  to  Ti  the  single  mask 

1  _  p-ie.2-‘^ 
t{u)  :=  ri(a;) - - - . 


Boxlets  correspond  to  the  choice  s  =  21  for  the  dilation  matrix.  The  cycles  of  21 
are,  of  course,  singletons,  hence  the  inductive  process  allows  us  to  insert  one  direction  per 
step.  We  may  begin  the  inductive  process  with  any  box  spline  (pi  whose  (integer)  direction 
set  is  a  basis  for  IR'^  (that  box  spline  is  then,  up  to  a  normalization  constant,  the  support 
function  of  some  parallelepiped).  As  the  initial  wavelet  masks  Ti,  we  take  the  tensor 
product  construction  (thus,  we  obtain,  up  to  a  linear  transformation,  the  multivariate 
Haar  wavelets).  We  select  any  sequence  B  =  2ind,  after  m  insertions  obtain 

the  final  box  spline 


${u})  =  $i{uj) 

i=i 


I  _ 

2i^j  •  u 
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There  are  2^^  -  1  +  m  boxlets  in  this  construction.  The  first  2^^  -  1  ones  correspond  to  the 
initial  tensor  product  masks  (which,  we  stress,  are  applied  to  the  smooth  0).  The  other 
masks,  which  we  index  by  H,  have  a  “triangular  structure”: 


J-1 

n 

i:=l 


1  +  e 


with  Ti  the  refinement  mask  of  the  initial 

Assuming  that  m  is  relatively  small,  most  of  the  boxlets  have  small  support.  For 
example,  if  we  wish  to  construct  boxlets,  then,  independently  of  the  spatial  dimension 
d,  we  may  do  with  m  =  3.  In  this  case,  there  are  2*^-1  wavelets  with  “small”  support, 
one  with  “large”  support,  and  two  wavelets  with  intermediate  size  of  supports.  □ 

We  remark  that,  in  the  above  construction,  we  may  start  with  any  refinable  (box 
spline)  (/.>i,  provided  that  a  derivation  of  a  tight  frame  from  <f)i  is  available. 

Example.  We  show  that  the  construction  of  §3  is  a  special  case  of  the  construction  of 
this  section.  Here,  is  a  four-direction  box  spline  (as  defined  in  (2.1))  with  multiplicities 
(0,m2,0,m2),  and  (f)2  is  a  four-direction  box  spline  with  multiplicities  (mi,0,mi,0)  (so 
that  each  box  spline  has  only  two  active  directions).  The  function  02  is  refinable  with 
respect  to  the  dilation  matrix  s  employed  in  §3,  and  with  mask  (up  to  an  exponential 
factor)  cos”^^  {uji/2).  Since  this  mask  is  univariate,  one  can  use  our  univariate  construction 
from  [RS2]  to  obtain  the  T2  masks 

As  to  01,  the  box  spline  0i  is  refinable  with  mask  (up  to  an  exponential  factor)  cos”^^  (a;2/2), 
and  hence  the  same  univariate  construction  can  be  repeated  (as  indeed  we  did) .  Otherwise, 
we  may  assume  by  induction  that  Ti,  the  wavelet  masks  associated  with  0i,  are  already 
given. 

Note  that  the  above  discussion  shows  that  the  construction  of  §3  could  have  been 
made  gradual:  starting  with  four-direction  box  spline  with  multiplicities  (0,  m2, 0,  m2),  we 
may  have  appended  the  four-direction  box  spline  with  multiplicities  (mi,0,  mi,0)  step  by 
step:  The  number  of  wavelets  will  be  then  unchanged,  (mi  -t-  m2)  but  a  certain  saving  in 
the  size  of  the  larger  support  wavelets  can  be  achieved  in  this  way.  □ 

Univariate  cubic  splines.  In  the  final  example  of  this  section,  we  employ  the  gen¬ 
eral  algorithm  of  this  section  in  the  construction  of  a  univariate  compactly  supported  tight 
frame  generated  by  cubic  splines.  Prior  to  doing  that,  we  recall  that  [RS2]  already  pro¬ 
vides  such  a  construction.  Its  construction  yields  four  generators  (^ji)^_i,  each  supported 
in  the  interval  [—2, 2]  whose  Fourier  transforms  satisfy 


0i(2w) 


cos^~^  (a;/2)  sin^'^^  (a;/2) 

0V2)^ 
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Figure  4.1.  The  four  piecewise-cubic  wavelets  that  were  con¬ 
structed  in  [RS2].  All  of  them  are  supported  in 
[—2, 2]. 


The  graphs  of  these  three  functions  are  given  in  Figure  4.1. 

Here,  we  construct  another  piecewise-cubic  tight  frame,  generated  by  3  wavelets, 
albeit  with  larger  support.  In  this  construction,  we  choose  as  the  Haar  wavelet,  and  02 
as  the  quadratic  B-spline.  Thus: 


01(W)  =  - ^ = 

ILO 

and  the  corresponding  masks  are 


Tl{oj)  = 


1  +  e' 


r2  =  rf . 


We  now  choose  Ti  to  be  the  singleton 


Ti  :={ 


1  —  e' 


2i 


(that  corresponds  to  the  Haar  wavelet),  and  choose  T2  to  consist  of  two  masks: 


T2  :={ 


1  —  e 
¥ 


\/3(l 


4i 


-}• 


It  is  obvious  that  Ti  satisfies  (b)  of  (4.1).  It  is  less  obvious,  but  still  can  be  checked 
directly,  that  T2  satisfies  (c)  of  (4.1).  Thus,  by  our  algorithm,  the  following  three  mother 
wavelets  generate  a  tight  frame  for  L2(1R): 

0i(t^)  =  - 


and 


,  1-f  e--/2x/3(l-e-2-)^^ 

=  — - - 

M^)  = - 2 - (  — 2i —  j 
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H — F 


Figure  4.2.  The  three  piecewise-cubic  wavelets  that  were  con¬ 
structed  here. 


where 


${u})  = 


The  supports  of  these  wavelets  are  [0,4],  [0,6],  [0,7],  respectively,  and  their  graphs 
are  shown  in  Figure  4.2.  We  note  that  "02  and  03  are  splines  with  integer  knots;  i.e.,  in 
standard  wavelet  terminology,  they  belong  to  Vq. 
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